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It is well known that waves propagating in a nontrivial medium develop "tails". However, the 
exact form of the late-time tail has so far been determined only for a narrow class of models. We 
present a systematic analysis of the tail phenomenon for waves propagating under the influence of 
a general scattering potential V{x). It is shown that, generically, the late-time tail is determined by 
spatial derivatives of the potential. The central role played by derivatives of the scattering potential 
appears not to be widely recognized. The analytical results are confirmed by numerical calculations. 
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One of the most remarkable features of wave dynam- 
ics in curved spacetimes is the development of "tails". 
Gravitational waves (or other fields) propagate not only 
along light cones, but also spread inside them. This im- 
plies that at late times waves do not cut off sharply but 
rather die off in tails. In particular, it is well established 
that the late-time evolution of massless fields propagat- 
ing in black-hole spacetimes is dominated by an inverse 
power-law behaviour. 

Price 0] was the first to provide a detailed analy- 
sis of the mechanism by which the spacetime outside 
of a (nearly spherical) collapsing star divests itself of 
all radiative multipole moments, and leaves behind a 
Schwarzschild black hole. The evolution of such waves 
(gravitational, electromagnetic, and scalar) in a curved 
spacetime is governed by the Klein-Gordon (KG) equa- 
tion @ 



^f' = 



(1) 



where the term 1(1 -I- l)/a;^ is the well-known centrifugal 
barrier (Z is the multipole order of the field), and Vix) 
is an effective curvature potential (we henceforth take 
Xs = 1 without loss of generality). It was demonstrated 
that all radiative perturbations decay asymptotically as 
an inverse power of time, the power indices equal 11 ~\- 3. 
Physically, these inverse power-law tails are associated 
with the backscattering of waves off the effective curva- 
ture potential, V^x), at asymptotically far regions [^[01. 

The analysis of Price has been extended by many au- 
thors. Bicak Q generalized the results to the case of an 
electrically neutral scalar field propagating in a charged 
Reissner-Nordstrom spacetime. Leaver demonstrated 
that the late-time tail can be associated with the ex- 
istence of a branch cut in the Green's function for the 
wave propagation problem. 

Gundlach, Price, and PuUin joj showed that power- 
law tails are a genuine feature of gravitational collapse 
- the existence of these tails was demonstrated in full 
non-linear numerical simulations of the spherically sym- 
metric collapse of a self-gravitating scalar field (this was 
later reproduced in 0). Moreover, since the late-time 
tail is a direct consequence of the scattering of the waves 



at asymptotically far regions, it has been suggested that 
power-law tails would develop independently of the ex- 
istence of an horizon j|] . This implies that tails should 
appear in perturbations of stars as well. Moreover, the 
KG wave equation has a much wider range of physical 
applicability. For instance, electromagnetic waves prop- 
agating in an optical cavity can be described by this 
equation as well |^,|lO| [A nontrivial dielectric constant 
distribution, (x) , provides a nontrivial medium, and is 
analogous to the presence of a scattering potential ^^(2^).] 
For other related works, see e.g., |l0|-|20||, and references 
therein. 

The dynamics of waves in rotating black-hole space- 
times has received much attention recently |2^|3^ . It has 
been demonstrated that power-law tails present in the 
Kerr spacetime as well (However, the damping indices are 
generally different from those found in the Schwarzschild 
case.) 

Yet, a thorough understanding of the fascinating phe- 
nomenon of wave tails is far from being complete. In 
particular, the exact form of the late-time tail has so far 
been determined only for a narrow class of models. In 
fact, most of previous analyses are restricted to "loga- 
rithmic potentials" of the form V{x) ~ In'^x/a;" (where 
a > 2 and /J = 0, 1 are parameters). In a brilliant work, 
Ching et. al. jl^ have shown that, generically, the late- 
time behaviour of waves propagating under the influence 
of these specific potentials has the form 'J ^ In'' 
(with the well-known exception of a "pure" power-law 
decay in the Schwarzschild spacetime). 

The purpose of this paper is to present a systematic 
analysis of the tail phenomenon with a general scattering 
potential V{x). Ching et. al. ||l^ provided a heuristic 
picture of the scattering problem which is an important 
first step in this direction: Consider a wave from a source 
point y. The late-time tail observed at a fixed spatial 
location, x, is a consequence of the wave first propagating 
to a distant point x' » y,x, being scattered by V{x'), 
and then returning to x at a time t :^ {x' — y) + {x' — a;) ~ 
2a;'. Hence, according to this scenario, the scattering 
amplitude (and thus the late-time tail itself) are expected 
to be proportional to V{x') ~ V{t/2). 

However, it is well-known, at least in the specific case of 
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logarithmic potentials, that this simple picture requires 
two modifications. First, there is an extra suppression 
of the late-time tail by a factor of Second, if a is 

an odd integer less than 2Z + 3, the leading term in the 
late-time tail vanishes , and one should consider sub- 
leading terms (The well-known Schwarzschild spacetime 
belongs to this case.) 

There are several interesting and important open ques- 
tions regarding the tail phenomenon in a general analy- 
sis. What determines the late-time tail - is it simply 
the asymptotic form of the scattering potential itself, as 
suggested by the heuristic picture ? How generic is the 
suppression of waves with ^ > 1 ? Is it always by a factor 
of ? What is the (most) general class of scattering 
potentials for which the leading term in the late-time tail 
vanishes ? (The Schwarzschild spacetime is only one spe- 
cific example to these potentials.) What is the general 
form of the late-time behaviour in these cases ? These 
questions, and several others call for a study of the gen- 
eral properties of wave tails. In this paper we present our 
main results for this fascinating phenomena. 

We consider the evolution of a wave field whose dynam- 
ics is governed by a KG-type equation -f V{x)^ = 
0. Resolving the field into spherical harmonics $ = 
'^\n{t T {6 T f) / r [r being the circumferential ra- 

dius), one obtains a wave equation of the form Eq. (|^) 
for each multipole moment | |3^ (For brevity we hence- 
forth suppress the indices l,m on ^'.) 

It proofs useful to introduce the double-null coordi- 
nates u = t — X and v = t + x, which are a retarded 
time coordinate and an advanced time coordinate, re- 
spectively. The initial data is in the form of some com- 
pact outgoing pulse in the range uq < u < ui, specified 
on an ingoing null surface v = vq. 

The general solution to the wave-equation (Q) can be 
written as a series depending on two arbitrary functions 
F and G |] 



k=0 
oo 

+ Y,Bi{x)\G^^-''-^'>iu) + {-lfF^^-''~^\v)\ , (2) 

fc=0 

wheieA[ = {l+k)\/2^k\{l-k)\. For any function i7, i/W 
is its fcth derivative; negative-order derivatives are to be 
interpreted as integrals. The functions B\.{x) satisfy the 
recursion relation 



B',^\[B'U-l{l + l) 



'-Bk-i-V{AkX-'' + Bk-i)\ , 

(3) 



for fc > 1, where B' = dB/dx, and B'^ = -V{x)/2. 

For the first Born approximation to be valid the scat- 
tering potential V{x) should approach zero faster than 



l/x^ as a; — > oo, see e.g., | |To| , [l6t . Otherwise, the scat- 
tering potential cannot be neglected at asymptotically 
far regions [see Eq. (|^) below]. It is useful to classify 
the scattering potentials into three groups, according to 
their asymptotic behaviour: 

• Group I: \V'\ approaches zero slower than |X^|/a;, 
and faster than |y| as a; ^ oo. 

• Group II: \V'\ approaches zero at the same rate as 
\V\ as X oo. 

• Group III: \V'\ approaches zero at the same rate as 
\V\/x as a; ^ oo. 

Group I. — For this case the recursion relation, Eq. 
(I), yields Bk{x) = -2-('=+i)T/('=-i) (cc). 

The first stage of the evolution is the scattering of the 
field in the region mq < u < ui. The first sum in Eq. 
(^) represents the primary waves in the wave front (i.e., 
the zeroth-order solution, with V = 0), while the second 
sum represents backscattered waves. The interpretation 
of these integral terms as backscatter comes from the 
fact that they depend on data spread out over a section 
of the past light cone, while outgoing waves depend only 
on data at a fixed u ||l| . 

After the passage of the primary waves there is no 
outgoing radiation for u > ui, aside from backscattered 
waves. This means that G{ui) — 0. Hence, for a large x 
a,t u — ui, the dominant term in Eq. (0) is 



*(m ui,x) ^ Bi{x)G^-^Hui) 



(4) 



This is the dominant backscatter of the primary waves. 

With this specification of characteristic data on u = 
Ml, we shall next consider the asymptotic evolution of 
the field. We confine our attention to the region u > ui, 
X ^ Xs- To a first Born approximation, the spacetime 
in this region is approximated as flat 0,^]. Thus, to first 
order in V (that is, in a first Born approximation) the 
solution for can be written as 



k=0 



(5) 



Comparing Eq. (|^) with the initial data on u = ui, Eq. 
d), one finds f{v) = -2-iy(-i) (w /2)G(-i) (ui). 

For late times t ^ x one can expand g{u) = 

oo 

J2 i-l)" g^"Ht)x"-/n\ and similarly for f{v). With these 

n=0 

expansions, Eq. (^) can be rewritten as 

oo 
n— — l 

where the coefficients if" are those given in They 
vanish for —l<n<lifl — nis odd. 
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Using the boundary conditions for small r (regular- 
ity as X — > — oo, at the horizon of a black hole, or at 
X = 0, for a stellar model), one finds that at late times 
g{t) = (— to first order in the scattering po- 
tential V (see e.g., for additional details). That is, 
the incoming and outgoing parts of the tail are equal in 
magnitude at late-times. This almost total reflection of 
the ingoing waves at small r can easily be understood on 
physical grounds - it simply manifests the impenetrabil- 
ity of the barrier to low-frequency waves (which are 
the ones to dominate the late-time evolution [^). We 
therefore find that the late-time behaviour of the field at 
a fixed radius {x <C t) is dominated by [see Eq. (^)] 



(2i + l) 



which implies 



(7) 



(8) 



where ^Pq = -2~(2'+i)is:^'+i. Hence, the late-time tail is 
determined by the 2fth derivative of the scattering po- 
tential. 

The analysis for Groups II and III proceeds along the 
same lines. That is, one should first solve Eq. (^ for 
Bk{x) find f{v) using Eq. (||) which finally yields 
'^{t 00 ) through Eq. (^). In the following we present 
the main results for groups I and II. 

Group II. — The dominant backscatter of the primary 

oo 

waves is ^'(w — ui,x) — ^ Bk{x)G'^^^^^'^\ui), where 

k=l 

the Sfe are the same as for group I. Using an analysis 
along the same lines as before, one finds 



n=i+l,i+3,... k=l 

xG('"'^-i)(Mi)V^("+'=-i)(t/2) , 



(9) 



at late-times. Note that Eq. (||) is merely a general- 
ization of Eq. (||), and reduces to it if \V'\ approaches 
zero faster than \V\ (in which case V^^^'^ dominates at 
late-times). 

Group III. — Let V{x) = W{x)/x"', where a > 2, 

and W{x)lx'^ ^ as x oo for any 7 > 0. The 

general solution for Bk{x) is now given by Bk{x) ~ 
00 

J2 b„{k,a)W^"Hx)/x°'+''-"-^, where the &„(fc,a) are 

71=0 

complicated numerical coefficients, which are determined 
by the recursion relation Eq. (^ . The analysis now pro- 
ceeds along the same lines as before; One finds that the 
late-time behaviour of the wave is given by 



* ~ G(~i)(wi)a;'+i ^c„(/,a)VK(")(i/2)t"(" 



-{a+2l-n) 



(10) 



TABLE I. Late-time behaviour for various scattering po- 
tentials. 



Group V{x) 



I < /3 < 1 

I sm{x'^)/x°', < /3 < 1 

II sin(a:)/a;" 



00) 



-(t/2)P 



-2i(l-/3) 



sin[(t/2)'']t-"-2i(i- 
periodic fun. xt^" 



Illa x^/" /3 > 

Illb a^i/^'^/s", a odd < 2; -I- 3, /? > t-(«+2'+'5) Int 



where the coefficients c„(Z,a) are constructed from the 
6„(/, a). 

Group III is divided into three subgroups according to 
the asymptotic behaviour of the function W{x), and the 
value of a: 

Subgroup Ilia: \W'\ approaches zero faster than | ly | /x 
as X — > 00, and a is not an odd integer less than 
21 + 3. In this case, the dominant term in Eq. ( [lO| ) 
is W{t/2)t-^"+^^\ RecaU now that V 
t — > 00, which implies 



W/t 



a+l 



as 



^ - G(-i)(ui)x'+il/(20(t/2) ^ 



(11) 



at late-times. 

Subgroup Illb: |II^'| approaches zero /asier than jH^j/x 
as X — > 00, and a is an odd integer less than 21 + 
3. (This subgroup of scattering potentials includes the 
Schwarzschild spacetime as a special case.) In this case 
one finds that the leading term of Bi (x) [proportional to 
V^^^^\x)] vanishes, and sub-leading terms should there- 
fore be considered. Hence, the late-time behaviour of the 
wave is dominated by 



*-G(-1)(mi); 



='+Vf)(t/2) 



(12) 



Here V^^''"' is the (first) sub-leading term in the asymp- 



totic derivative. Namely, vjf'^ ~ W /t°'+'^'^-'^ . [Note 
that the results of |l^ for the specific family of loga- 
rithmic potentials (with W ~ In'^x), coincide with the 
general expressions, Eqs. (|ll] ) and (p^.j 

Subgroup IIIc: \W'\ approaches zero at the same rate 
as |iy|/x as X — * 00. In this case, the late-time dynamics 
of the field is given by Eq. ([l0|). 

Numerical calculations. — It is straightforward to in- 
tegrate Eq. (|l|) using the methods described in [p|,p^. 
The late-time evolution of the field is independent of the 
form of the initial data used. The results presented here 
are for a Gaussian pulse. 

Table | gives a selected list of scattering potentials, 
chosen as representative of the various different groups 
(We have studied other potentials as well, which arc not 
shown here.) 

The temporal evolutions of the waves (under the in- 
fluence of the various scattering potentials) are shown in 
Figs. and ^. We find an excellent agreement between 
the analytical results and the numerical calculations. 
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FIG. 1. Temporal evolution of the field for scatter- 
ing potentials which belong to Group I. Top panel: 
V{x) = e~'"^/'^°' (the results presented here are for the par 
rameters (5 = 1/2 and xo = 1/2). We display the quantity 
ln|^'t2'(i-'5)| vs t". There is a definite linear dependence at 
late-times, with slopes of —0.99 and —1.02 for I = (lower 
graph) and / = 1, respectively. These are in excellent agree- 
ment with the analytically predicted value of —1. Bottom 
panel: V{x) = sm[{x / xq)^] / x°' (the results presented here 
are for the parameters a = 3, /3 = 1/2, and xo = 1/2). An 
oscillatory power-law fall off is manifest at late times. The 
power-law indices (determined from the maxima of the oscil- 
lations) are —4.05, and —5.07, for Z = 1 (upper graph), and 
1 = 2, respectively. These are in excellent agreement with the 
analytically predicted values of —4, and —5. 

Summary and physical implications. — We have given 
a systematic analysis of the tail phenomena for waves 
propagating under the influence of a general scattering 
potential. 

It was shown that the late-time tail is governed by spa- 
tial derivatives of the scattering potential (generically, by 
the 2Zth derivative). In particular, the potential function 
itself does not enter into the expression of the late-time 
tail (with the exception of the monopole case). The cen- 
tral role played by derivatives of the scattering potential 
appears not to be widely recognized. The analytical re- 
sults are in excellent agreement with numerical calcula- 
tion. 

In addition, wc have demonstrated that the (extra) 
suppression of waves by a factor of (which adds to 
the basic late-time decay), a phenomena well-known in 
black-hole spacetimes, is actually not a generic feature 
of the scattering problem. In particular, for scattering 
potentials that belong to group I the suppression of the 
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FIG. 2. Top panel: Temporal evolution of the field for a 
scattering potential of the form V{x) = sin(a;/a;o)/a;", which 
belongs to Group II (the results presented here are for the 
parameters a = 3 and xo = IOO/tt). The power-law indices 
are —3.02 and —3.04 for / = 1 (upper graph) and I = 2, re- 
spectively. These values are in excellent agreement with the 
analytically predicted value of —3 for both I = 1 and I = 2. The 
oscillations period agrees with the predicted value to within 
0.4%. Bottom panel: Temporal evolution of the field for scat- 
tering potentials of the form V{x) = x^^'^ /x", which be- 
long to Group III (the results presented here are for the case 
(5 = 1). We present the quantity In j^t" In'' t\ vs t, where: (i) 
e = 3, (5 = for a = 3 with I = 0, (ii) e = 6, 5 = — 1 for a = 3 
with 1 = 1, and (iii) e = 5 = Q for a = A with I = 1 (from 
bottom to top). These quantities are analytically predicted 
to approach a constant value at asymptotically late-times. 

waves is smaller, while for scattering potentials that be- 
long to group II there is no (extra) suppression at all. 

Moreover, it was shown that the familiar case of the 
Schwarzschild spacetime belongs to a wider group of 
scattering potentials, in which the leading term in the 
tail [proportional to V^'^''\t/2)] vanishes (and thus, sub- 
leading terms dominate the late-time dynamics) . 

We are at present extending the analysis to include: (i) 
time-dependent scattering potentials, and (ii) scattering 
potentials that lack spherical symmetry (in which case 
the scattering problem is of 2 -|- 1 dimensions). 
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